Zero-energy states of graphene triangular quantum dots in a magnetic field 
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We present a tight-binding theory of triangular graphene quantum dots (TGQD) with zigzag edge and broken 
sublattice symmetry in external magnetic field. The lateral size quantization opens an energy gap and broken 
sublattice symmetry results in a shell of degenerate states at the Fermi level. We derive a semi-analytical form for 
zero-energy states in a magnetic field and show that the shell remains degenerate in a magnetic field, in analogy 
to the 0th Landau level of bulk graphene. The magnetic field closes the energy gap and leads to the crossing of 
valence and conduction states with the zero-energy states, modulating the degeneracy of the shell. The closing of 
the gap with increasing magnetic field is present in all graphene quantum dot structures investigated irrespective 
of shape and edge termination. 
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I. INTRODUCTION 



Graphene currently attracts considerable attention due to 
remarkable electronic and mechanical properties. 1-11 When 
graphene is reduced to graphene nanostructures, new effects 
related to size-quantization and edges appearJi^ Consider- 
able experimental effort has been made aiming at producing 
graphene nanostructures with desired shape and edges^^— 
Among graphene nanostructures, nanoribbons and quantum 
dots are of particular interest. In graphene quantum dots, a 
size-dependent energy gap opens,— ^ and its magnitude is 
determined by shape and edge termination. In graphene quan- 
tum dots with zigzag type edges edge states with energy in 
the vicinity of the Fermi energy appear . 5 i 33 ' 36 ^ These edge 
states have significant effects on low-energy electronic prop- 
erties such as a decrease of the energy gap compared to struc- 
tures with armchair termination or, when combined with bro- 
ken sublattice symmetry, a creation of the degenerate shell of 
zero-energy states in the middle of the energy gap*2MM2 

The influence of an external magnetic field on the elec- 
tronic properties of the graphene quantum dots was also 
studied ] 18 ' 34,48 ^ The magnetic field plays the role of a 
tunable external parameter allowing to change electronic 
properties in a controllable way. Graphene quantum dots 
and rings with circular, square, hexagonal, triangular, and 
rhombus-shaped shapes with zigzag and armchair edges were 
investigated ) ' 50 ~ 52 ' 54 ^ Triangular graphene quantum dots 
with reconstructed edges, consisting of a succession of pen- 
tagons and heptagons, were also considered. 59 The com- 
parison between tight-binding and continuum model, the 
Dirac-Weyl equation, was analyzed for graphene quantum 
dots with different type of edges: zigzag, armchair, and 
infinite-mass boundary conditions. — — For a circular dot, 
good qualitative agreement between experiment and ana- 
lytical model with infinite-mass boundary conditions was 
obtained^^ 9 . Magneto-optical properties were also theoreti- 
cally investigated. 34 ' 56 The absorption spectra differ for hexag- 
onal structures with armchair and zigzag edges due to differ- 
ent level structures and the oscillator strengths. A fast reduc- 
tion of the energy gap with increasing magnetic field in zigzag 
hexagon in comparison with zigzag triangle was noted^ 6 ^ 



In this work, we present a tight-binding theory of triangular 
graphene quantum dots(TGQD) with zigzag edge and broken 
sublattice symmetry in external magnetic field. The lateral 
size quantization opens an energy gap and broken sublattice 
symmetry results in a shell of degenerate states at the Fermi 
level. Building on our previous work 44 we derive here a semi- 
analytical form for zero-energy states in a magnetic field and 
show that the shell remains degenerate at all magnetic fields 
perpendicular to the plane of the TGQD, in analogy to the 
0th Landau level of bulk graphene. However, we find that the 
magnetic field closes the energy gap and leads to the crossing 
of valence and conduction states with the zero energy states, 
modulating the degeneracy of the shell. The closing of the 
gap with increasing magnetic field is present in all graphene 
quantum dot structures investigated irrespective of shape and 
edge termination. 

The paper is organized as follows. In Sec. II, we present a 
brief outline of the tight-binding model with an incorporation 
of a perpendicular magnetic field. The analysis of the evolu- 
tion of the energy spectra of TGQD, a derivation of the ana- 
lytical form for eigenfunctions corresponding to zero-energy 
states, and a prediction of crossings of valence and conduc- 
tion states with E = are included in Sec. III. In Sec. IV 
the energy gap in a magnetic field for GQDs with different 
shapes and edge termination is considered. The conclusions 
are presented in Sec. V. 



H. MODEL 

We describe graphene quantum dots using the nearest- 
neighbor tight-binding model which has been successfully 
used to describe graphene 1 and applied to other graphene 
materials such as nanotubes, nanoribbons and quantum 
dot a 33 ' 36,39 ^ 2,44 ' 60 . A perpendicular magnetic field can be in- 
corporated by using Peierls substitution 61 . The Hamiltonian 
reads, 
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FIG. 1 : Left: Energy spectrum of triangular graphene quantum dot 
with N = 97 atoms and N e d = 7 degenerate zero-energy states in 
the absence of a magnetic field. Right: the evolution of the spectrum 
from the left in a magnetic flux. The degenerate zero-energy shell is 
immune to the magnetic field. The states from the conduction and 
valence bands, labeled by I and II, respectively, crosses for (p/<po — 
0. 1 1, closing the energy gap. 



where t is hopping integral, a! (b. ) and a !£J -(&;o-) are creation 
and annihilation operators on a site i corresponding to sublat- 
tice A(B) of bipartite honeycomb lattice, (2, j) indicate sum- 
mation over nearest-neighbors, and cr is spin index. Hopping 
integral between nearest neighbors is t — -2.8 eV 11 . Under 
symmetric gauge, a vector potential A = B z /2(-y, x, 0), and 



2n — I Ac/1 = 2n — — fxy,- — JC,vA (2) 
he J,, 2cf>o V ,yj in > 



corresponds to a phase accumulated by electron going from 
site i to j, which is equal to a magnetic flux going through 



area S = 



spanned by vectors r, and rj, and 4>a — — is 



magnetic flux quantum. The evolution of the energy spectrum 
in a magnetic field will be shown in units of the magnetic flux 
threading one benzene ring, <f>/</>o = B Z S o/<Po, where So = 
3 V3fly/2 is benzene ring area with ao - 1.42 A. 



III. ZIGZAG TRIANGULAR QUANTUM DOT IN A 
MAGNETIC FIELD 

A. The evolution of the energy spectrum 

We focus here on the effect of the magnetic field on the 
electronic properties of TGQDs, quantum dots with broken 
sublattice symmetry. We illustrate the energy spectrum and 
its evolution with increasing magnetic field on a TGQD with 
N - 91 carbon atoms. Fig. Q] shows the energy spectrum 
and it's evolution in the magnetic field obtained by numerical 
diagonalization of the Hamiltonian, Eq. {T). 

At B — there are Ndeg = 7 degenerate states at zero energy 
or Fermi level. The number of states is equal to the difference 
between the number of A and B atoms 44 . The states belonging 
to the degenerate shell are primarily localized at the edge of 
the triangle and are entirely localized on one sublattice, say A, 
as shown in Fig. |2c). 



The evolution of the energy spectrum as a function of the 
magnetic field is shown on the right hand side of Fig. Q] 
The spectrum is symmetric with respect to E — due to 
electron-hole symmetry. This symmetry is broken when hop- 
pings to the second nearest neighbors in Hamiltonian, ((T), are 
included. The highest valence state and the lowest conduction 
state with E = +1.57 eV, which in the absence of the mag- 
netic field are each doubly degenerate, split in the presence of 
a magnetic field. The state labeled by II from the valence band 
increases and the state labeled by I from the conduction band 
decreases its energy with increasing magnetic field, closing 
the energy gap. Around (f>/<f>o - 0.1 1 these states reach Fermi 
level at E = 0. 

The explanation of why the energy gap closes in a mag- 
netic field can be found by considering Dirac Fermions in bulk 
graphene^i. We focus on one of two Fermi points, say K 
point. Following Refs. 62 63 the energy spectrum of Dirac 
Hamiltonian in the presence of magnetic field is given by 



E = ±^2Hv F eB z \n\/c, 



(3) 



where v F is Fermi velocity, c speed of light, and n Landau 
level index. The ± sign corresponds to electron (hole) Landau 
levels. A unique property of the energy spectrum is the exis- 
tence of the n — Landau level with energy E — 0, constant 
for all magnetic fields. When the magnetic field is applied to 
graphene quantum dots, discrete energy levels evolve into the 
degenerate Landau Levels for Dirac Fermions. Thus, some 
levels have to evolve into the 0-th Landau level, closing the 
energy gap as shown in Fig. Q] Another feature of the 0-th 
Landau level is that the wavefunctions are localized on only 
one sublattice, similarly to the zero-energy states in TGQD 44 . 

We note in Fig. [T]that the zero-energy degenerate shell is 
immune to the magnetic field as is the n — Landau level. 
This is certainly different from electronic states in semicon- 
ductor quantum dots, where ~ B 2 dependence is observed. 64 

These comments are now illustrated by examining wave- 
functions of a TGQD in a magnetic field. We investigate the 
evolution of the probability density of the wavefunction corre- 
sponding to state I, bottom of the conduction band, from Fig. 
[1] and the total probability density of the zero-energy degen- 
erate shell in a magnetic field. For state I, probability densities 
at low and high magnetic field values are shown in Fig. |2ja) 
and (b), respectively. We note that due to the electron-hole 
symmetry, an identical evolution for the state II from Fig. [1] 
(not shown here) occurs. Eigenfunctions of states with energy 
-\E\ and +\E\ differ only by a sign of a coefficient on sublat- 
tice B indicated by filled circles in Fig. [2] giving identical 
electronic densities. For <pl<po — 0.01, Fig. |2ja), the state 
I is mostly localized at the center of the dot. With increas- 
ing magnetic field, it starts to occupy the edge sites, shown 
for <pl4>Q — 0.08. We note that for arbitrary magnetic field 
this state is equally shared over two sublattices, i.e, has 50% 
sublattice content. In Fig. |2|c) and (d) the evolution of the 
total electronic density of the degenerate zero-energy shell is 
shown. The electronic density of the degenerate shell is ob- 
tained by summing over all Ndeg - 7 states. Initially, degener- 
ate states are strongly localized on edges, shown in Fig. |3Jc) 
for <pl<po — 0.01. When the magnetic field increases, these 
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FIG. 2: (Color online) The evolution of electronic densities in a mag- 
netic field of the lowest state from the conduction band labeled by /, 
(a) and (b), and the degenerate shell of Ndeg = 7 energy levels (the 
sum of electronic densities of all Ndeg = 7 degenerate states), (c) and 
(d). The radius of gray circles is proportional to the electronic proba- 
bility density on a given site, (a) The state / for <p/(po - 0.01 is mostly 
localized in the center of the dot and with increasing magnetic field 
starts to occupy also edges, shown for (/>/t/>o - 0.08 in (b). (c) The 
degenerate states for (p/ipo — 0.01 are strongly localized on edges and 
for 4>I4>q - 0.08 move slightly to the center of the triangle, shown in 
(d). 

states move slightly towards the center of the triangle, shown 
in Fig. [2jd) for <f>/<f>o - 0.08. We note that even in the pres- 
ence of an external magnetic field states from the degenerate 
shell are still localized on only one type of atoms, sublattice 
A, indicated by open circles in Fig. [2] 

B. Analytical solution for zero-energy states 

Fig. Q] shows that numerical diagonalization of the tb- 
Hamiltonian gives the zero-energy states immune to external 
magnetic field. We will now prove this analytically. Our first 
goal is to show the existence of and find an expression for 
zero-energy eigenstates in the presence of a magnetic field. 
The zero energy states, if they exist, must be solutions of the 
singular eigenvalue problem, 

H TB *¥ = 0, (4) 

where the Hamiltonian Hjb is given by Eq. (Q~|). There is 
no coupling between two sublattices and the solution can be 
written separately for A-type and B-type of atoms. We first 
focus on sublattice A with an eigenfunction given by 

|T A ) = J] Cja)\0) = J] Cj\<pj), (5) 

j j 

where C, are expansion coefficients of eigenstates written in 
a basis of p z orbitals 04 localized on A-type site j for either 
spin state ommitted in what follows. 



According to Eq. |4] the coefficients Cj corresponding to 
one type orbitals localized around the second type site i obey 

where the summation is over j-th nearest neighbors of an atom 
i. In other words, the sum of coefficients multiplied by a phase 
e "Pji gained by going from one type site i to the other type site j 
around each site ; must vanish. For the z'-th B-type site plotted 
on the left in Fig. [3] Eq. <JSJ gives 

Cje^J + C k e^' k + Cie mx = 0, (7) 

where phases qjy, cp*, tfn are given by Eq. (|2j. Using a fact 
that qpa = -(pki for arbitrary i and k, Eq. (0 can be written as 

Ci = - (Cje-^ 1 + C k e- itftl ) , (8) 

where qpy; = cp ;7 + cp ;; and tp^ = tp*; + (fa correspond to phase 
changes going from A-type sites k to j, and I to /', respectively, 
through B-type site i (see right part in Fig. |3). Thus, in anal- 
ogy with the zero magnetic field case 44 , a coefficient from a 
given row can be expressed as a sum of two coefficients from 
an upper lying row, Cj and Q on the right in Fig. 3. The effect 
of the magnetic field is incorporated in the extra phase gained 
by going from a given site from an upper row of atoms to a 
lower one. For a reason which will become clear later, instead 
of using indices i, each A-type site will be labeled by two inte- 
ger numbers, i = {«, m). The first index, n, corresponds to an 
atom number in a given row counted from left to right, and the 
second one, m, corresponds to the row number. Let us illus- 
trate our methodology on a hexagonal benzene ring with three 
auxiliary A-type atoms with indices Co,o, £2,0 and Co.i, shown 
in Fig. Ufa). Eq. © can be used to obtain coefficients Cqj 
from Co,o and Ci,o, and Ci,i from C^o and C^o- Next, using 
Co,i and C\$, one obtains coefficient Cq,z, 

C ,2 = Co,oe- /cp ' + C,,o (e-^ + + Cy*" ** (9) 

with phase changes cp,, i = 1,2, 3,4, shown as black arrows in 
Fig. 13 a). The paths related to phase changes cp, go through in- 
termediate atomic sites, e.g., for cpi the path goes from a site 
Co to C,i through an intermediate fi-type atomic site, and 
next from a site C,i to C,2 through connecting fi-type atomic 
site. According to Eq. (0 and Fig. |4ja), there is one path con- 
necting Co and C,2> ° ne connecting Co an d ^0,2, but there 
are two paths around a hexagonal benzene ring connecting co- 
efficients Co an d Co,2- We have shown that the coefficient in 
the bottom, C2, can he expressed as a linear combination of 
coefficients from the top row, Co- ^ e w iU now demonstrate 
that all coefficient in arbitrary size triangles can be expressed 
in terms of coefficients Co- 
in Fig- Efb) a small triangle with N e d = 2 atoms on the 
one edge is plotted. Three auxiliary atoms with coefficients 
C,o> Co> an d C,3 were added. The total number of atoms is 
N = 16. In a similar way to the procedure used to obtain Eq. 
(0, a coefficient C,3 ca n be expressed as a sum of coefficients 
from the top. Here, from coefficients Co an d Co to C3 there 
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FIG. 3: (Color online) Left: phase changes going from S-type site i 
to three nearest neighbors A-type sites j, k, and /. The sum of three 
A-type coefficients multiplied by corresponding phase changes must 
vanish for zero-energy states. Right: phase changes going from A- 
type site k to /, and / to j, A coefficient from the bottom, C/, can 
be expressed as a sum of coefficient from an upper row, C, and Q, 
multiplied by a corresponding phase changes. 



(a) 



(b) 




FIG. 4: (Color online) (a) A hexagonal benzene ring with three aux- 
iliary comer atoms added. Each A-type site is described by two num- 
bers {n,m}. Black arrows indicate phase changes related to the paths 
going from an upper row of atoms,with indices fn, 0}, to an atom 
from in the bottom, Co,2- (b) Triangular zigzag graphene quantum dot 
with N ed = 2 atoms on the one edge. There is one path going from 
site Co.o to Co,3, and three paths going from site C\$ to C03. The num- 
ber of paths connecting a site {11, m) with a site from the top {n + j, 0) 
can be described by binomial coefficient N pat h{n,m,n + j) = 
< j <m. 



is only one path for each coefficient, and three paths for each 
coefficient connecting C^o to Cn,3, and C^o to C03. For trans- 
parency, only for the first two coefficients from the left (Co,o 
and C10) paths are plotted in Fig. Hfb). The number of paths 
from a given site in the upper row of atoms to lower lying 
atomic sites corresponds to numbers from a Pascal triangle, 
{1, 2, 1} for coefficient Cq,2, shown in Fig. Ufa), and {1,3, 3, 1} 
for coefficient Cn,3, shown for the first two coefficients from 
the left in Fig. Etb). The number of paths connecting a site 
{«, m] with a site from the top \n + j, 0} can be described by 
binomial coefficient N pat h(n,m, n + j) = < j < m. The 
general form for an arbitrary coefficient expressed in coeffi- 



cients from the top row can be written as 

m (./) 

C„^ = (-) m 2Z c ^.o e "* +/( °' 

j=0 i=i 



(10) 



where two numbers n and m satisfy condition < n, m < 
N e d + 1, and qp n +j(0 is a path-dependent phase change from 
a site \n + j, 0} to {«, m). One can note that in the absence of 
a magnetic field cp n +X0 = and Eq. ( TTOb reduces to Eq.(2) 
from Ref4i. 

The summation over all possible paths in Eq. ( TTOb is not 
practical. We now show a way of reducing the number of 
paths to only one. We use the fact that a phase change corre- 
sponding to a closed path around a hexagon is by definition 
qp c = 27r0/0o- The sum of two exponential terms standing 
next to coefficient Ci,o in Eq. (O can be written as 



-tm + e -m = (gito-q*) + i} e -m = + A e -^ t 



(11) 



where qp3 - cp2 = 2n(/)/<po is a closed path around a single 
hexagon, see Fig. Eta). Similarly for three exponential terms 
corresponding to paths connecting C10 and C03, shown in Fig. 
Hfb), one can write 



-(■ 



e^Pfr + e^ti (12) 



where (qp 4 - - qp 2 <) circles two hexagons and (cp4< - (P3O only 
one, see Fig. |4fb). Note that phases qp3 in Eq. (fTTT i and qj^ 
in Eq. (fT2l correspond to paths going on the right edge of the 

triangle. The sum of exponential terms of type e 2m(,/ *o ) with j- 
integer in Eq. (fTTT i and Eq. (fT2l forms geometric series which 
can be written as 



k 

Z- 

7=0 



J *o = 



1 



1 - e *o 



(13) 



with k determined by the number of encircled benzene rings, 
and k + 1 = (j \ is a number of paths connecting site {« + j, 0} 
to {n,m}, k = 1 in Eq. ( fTTT i and k — 2 in Eq. 0121 . see Fig. 
|U Using Eq. (fT3l . the number of paths in Eq. (TTOb can be 
reduced to only one. Eq. (TTOb can be written as 



m , 2^!('")t- 

^'n+ 7,0 ' 

7=0 



, 2ni ±- 

l - e ^ 



(14) 



where cp„ + j is the phase corresponding to the path on the right 
edge connecting site \n + j,0} and \n,m}. The coefficients 
C n ,m for all A-type atoms in the triangle are expressed as a lin- 
ear combination of coefficients corresponding to atoms on one 
edge, i.e., C,;o- There are N e d + 2 coefficients in an upper row 
of atoms, Cjfi, with < j < N e d + 1, which gives N e d + 2 inde- 
pendent solutions. Applying three boundary conditions corre- 
sponding to auxiliary atoms, C ,o = C/v^+i.o = Co^+i = 0, 
leaves only N e d - 1 solutions, which corresponds to the num- 
ber of zero-energy states, similar to the result obtained in the 
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(b) 




FIG. 5: (Color online) (a) A hexagonal benzene ring with coefficients 
d assigned to each B-type site. Black arrows indicate phase changes 
related to the paths going from one B-type site to another, (b) Trian- 
gular zigzag graphene quantum dot with N e j = 2 atoms on the one 
edge. Black arrows indicate phase changes related to the paths going 
from one B-type site to another along triangular three edges. White 
arrows indicate phase changes related to the paths going through the 
center. 



absence of a magnetic field in Ref. 44, We note that the so- 
lutions given by Eq. ( fT4l > are smooth functions of magnetic 
field, and exist for any value of (p. Thus they do not include 
zero-energy solutions corresponding to the crossing of con- 
duction and valence states with E — 0, e. g., for <f>/<f>o - 0.11 
for the triangular dot with N e d - 8 and N = 97 atoms, see Fig. 
Q] We investigate this issue by analyzing B-type atoms. 



C. Prediction of crossings of valence and conduction states 
with E = 

We consider the solution of Eq. <j4j> corresponding to wave- 
function localized only on B-type atoms. In Fig. [5] the same 
structures as in Fig. |4] without auxiliary corner atoms are 
shown with coefficients assigned to B-type atoms. For sim- 
plicity, only one index for each coefficient is used. According 
to Eq. [6] for a benzene ring plotted in Fig. Qa) we can write 



c 2 


= -C 1 e-" Pl2 , 


(15) 


c 3 




(16) 


C] 




(17) 



where phase changes from site i to j, qpy, are indicated in Fig. 
|5Ja). Eq. (15[ can be substituted into Eq. (fTrjj l. and next Eq. 
(fTBT l into Eq. (fTTT i. giving 



Ci =Ci(~i)V i(<Pi2+<te+<,>3l) , 



(18) 



which is satisfied for arbitrary C\ . Eq. [T8]leads to a following 
condition 



CP12 + <P23 + Cf31 + 7T = 2nk, 



with k = 0, +1,±2, .... A phase change in Eq. ([T9T l corre- 
sponds to a closed path around a single hexagon, tpi2 + 923 + 




FIG. 6: Energy spectrum as function of the magnetic flux for (a) a 
hexagonal benzene ring and (b) a triangular graphene quantum dot 
with Ned = 2 atoms on the one edge and a total of N = 13 atoms. A 
crossing of conduction and valence states with E = occurs for (a) 
<p/(po = 1/2 and (b) 0/0 o = 1/3 and (f>/(p = 2/3. 



cp3i = 2n(f>/(f>Q. A condition for crossing of a valence and con- 
duction states with E = is 



= 1/2. 



(20) 

In order to confirm validity of Eq. (l20t . we show the energy 
spectrum of a benzene ring as a function of a magnetic field 
in Fig. |6fa). The crossing of energy levels at E = occurs for 
4>l<f>o = 1/2, in agreement with Eq. (120) . 

We carry out a similar derivation for triangular zigzag 
graphene quantum dot with /V = 1 3 carbon atoms and Ned - 2 
atoms on the one edge, shown in Fig. EJb). A coefficient from 
the left upper corner, C\ , determines a coefficient Ci, 

C 2 = -de- i(f12 . (21) 

Next, a coefficient C3 can be determined by a coefficient Ci 

C 3 = -C 2 e~ itte (22) 

and combining with Eq. (f2T~b gives 

C 3 = (-l) 2 CV-' Xcpi2+tp23) = (~l) 2 Cie-* Pl3 . (23) 

Going in this way along the three edges of the triangle a closed 
loop, shown with black arrows in Fig. |3b), can be created. In 
the case of N e d = 2 shown in Fig. |3b), one goes through 
all B-type coefficients, while in larger triangles one goes only 
through outer coefficients. Thus, all outer B-type coefficients 
can be expressed by one chosen coefficient, C\ in this case. 
The loop from Fig. |5jb) can be written 



d = (-1) 6 CV 



-i'(cPl6+<P6l) 



= cv 



6m'-2m'(30/0o) 



(24) 



where the phase change qpie = 2; cp/,i+i> an d we use d a 
fact that the total phase change corresponds to a closed loop 
around three benzene rings, cpi6 + (f6i = 2n(3(p/(f>o)- Eq. d24l) 
gives a condition 



(19) and finally 



2kn = 6tt — 6n<pl4>o, 



3-k 



(25) 



(26) 



6 



with k = 0, ±1, ±2, .... Eq. d26b can be extended to different 
size triangles. The number of benzene rings in a triangle is 
Nb = N e d{Ned + l)/2, and Eq. (f26t can be written as 



0/00 



3JV erf - 2k 

Ned (Ned + 1)' 



(27) 



For the triangle with N e d = 2, Eq. d26j > predicts crossings for 
0/0o = 0, ±1/3, ±2/3, 1, ..but according to Fig. |6fb) there are 
no crossings for 0/0o = and 0/0o = 1. This is related to 
an extra condition in the center of the triangle, for coefficients 
C2, C\, and Ce- Phase changes between these coefficients are 
indicated by white arrows in Fig. |5}b). We can write 



and also 



C 6 = - (C 2 e- i(f26 + C 4 e- itfi6 ) 



C 6 = -Cie- i( - (p6l) = -Cie icp61 , 

c 2 = -c ie -^\ 

C 4 = (-ifCie-'f", 



(28) 



(29) 



where the phase change 914 = cp/,/+i ■ Combining Eq. d28| i 
and Eq. ( 1291 we get 

- CV'P 61 = - (-Cie- i(cpi2+cp26) + (-l) 3 Cie" ,Xcpl4+tp46) ) , (30) 

which gives 

_ l — g-*'(<Pt2+<P26+<P6l) + g-'(<Pl4+<P46+<P6l) (3 ^) 

With help ofFig.|5|b), we can notice cpi2 + cp26 + <P6i = 2n<P/<f>o 
and cpi4 + qp46 + cp6i = 2n(2(p/<f)Q). Thus, we can write 



1 + e -2ni<p/4>o + e -2m(2<p/4> ) _ q 



or using a sum of geometric series 



I _ g-2m(30/0o) 



- e 



(32) 



(33) 



Eq. (f33l > gives a solution for -2n(3(p / <po) = 2nk, £- integer, 
and finally 0/0o = —k/3, but with an extra condition 0/0o + I, 
with I = 0, ±1, ±2, .. due to a denominator. This is in agree- 
ment with Fig. |6|b). We note that for all triangles, the predic- 
tion of crossings of conduction and valence states with E = 
given by Eq. d27| > has to be supported by extra conditions from 
equations for coefficients from the center of the triangle. For 
example, for the triangle with N e d = 8, the first crossing oc- 
curs for <f>/4>o = 1 19, while incomplete condition given by Eq. 
d27| i predicts the first crossing for <f>l<fro = 1 /36, and the fourth 
crossing for <f>/<f>o = 1 /9. 

This concludes the analysis of the magnetic field depen- 
dence of the energy spectrum of triangular graphene quantum 
dots with zigzag edges. 



IV. SHAPE AND EDGE DEPENDENCE OF THE ENERGY 
GAP IN A MAGNETIC FIELD 

We discussed above the magnetic field closing of the en- 
ergy gap in triangular graphene quantum dots. In Fig. [7] we 
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FIG. 7: The energy gap as function of the magnetic flux for triangu- 
lar zigzag quantum dot with N = 622 atoms (solid line), hexagonal 
zigzag quantum dot with N = 600 atoms (dotted line), hexagonal 
armchair quantum dot with N = 546 atoms (dashed line). 



analyze the evolution of the energy gap in graphene quan- 
tum dots with different shapes and edges in a perpendicular 
magnetic field. The energy gaps as a function of a mag- 
netic field obtained by diagonalizing Hamiltonian given by 
Eq. (Q]) are shown for three different types of quantum dots; 
zigzag triangle, zigzag hexagon, and armchair hexagon. All 
three structures have similar sizes, consisting of N 600 
atoms with area 5-14 nm 2 . The energy gap corresponds 
to the difference between the energy of the lowest state from 
the empty conduction states and the highest state from the 
doubly occupied valence states. In the absence of magnetic 
field, the zigzag triangular graphene quantum dot has a sig- 
nificantly larger gap then for hexagonal armchair and zigzag 
dots. For zigzag triangle and armchair hexagon, the energy 
gap decays as the inverse of the square root of the number 
of atoms E gap ~ V/Y 35 . When the magnetic field increases, 
the energy gap closes for all structures. The closing of the 
gap is related to the presence of the 0-th Landau level in bulk 
graphene. The bulk limit is reached in a very large magnetic 
field. The closing of the gap gives an opportunity to control 
optical properties of graphene quantum dots with the external 
magnetic field. 



V. CONCLUSIONS 

The electronic properties of triangular graphene quantum 
dots with zigzag edges and broken sublattice symmetry in 
the presence of perpendicular external magnetic field were 
described. It was shown that the degenerate shell of zero- 
energy states in the middle of the energy gap is immune to 
the magnetic field in analogy to the n — Landau level of 
bulk graphene. An analytical solution for zero-energy states 
in the magnetic field was derived. The energy gap was shown 
to close with increasing magnetic field, reaching zero at spe- 
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cial values of the magnetic field. The gap closing was found 
independent of quantum dot size, shape, and edge termination. 
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